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1 .  Introduction/Background 


The  Ballistic  Research  Laboratory-Computer- Aided  Design  (BRL-CAD)  is  a  solid  geometry 
modeling  CAD  system  that  assists  the  users  in  constructing  and  analyzing  complex  structures 
that  would  otherwise  require  advanced  computations.1,2  BRL-CAD  uses  a  tool  called  the  ray- 
tracer3  to  help  analyze  solid  objects.4,5  This  tool  has  the  ability  to  quickly  render  graphics  and 
perform  various  operations  pertinent  when  analyzing  such  geometries  as  BoTs  (Bag  of 
Triangles).6  A  BoT  is  generally  used  to  represent  triangle  mesh  objects  in  BRL-CAD.  BoT 
objects  have  several  modes,  such  as  surface  mode,  volume  mode,  and  plate  mode.  The  mode  that 
was  primarily  used  throughout  this  project  is  the  plate  mode,  which  is  used  to  represent  a  solid 
object  in  the  form  of  a  sheet  with  a  small  uniform  thickness. 

The  goal  of  this  project  was  to  create  all  possible  triangles  of  variegated  designs,  lengths,  and 
angles  within  a  specified  finite  range  and  then  to  perform  tests  to  assess  the  functionality  of  the 
ray -tracer  on  these  objects.  These  tests  were  used  to  gather  information  on  whether  the  ray-tracer 
can  hit  points  at  various  locations  relative  to  each  triangle,  such  as  hitting  directly  on  the  triangle 
or  just  slightly  exterior  to  the  triangle.  The  ability  to  extensively  test  the  BRL-CAD’s  ray-tracer 
on  BoTs  will  be  useful  in  identifying  the  baseline  accuracy  of  the  ray-tracer  and  then  in 
identifying  ways  to  improve  its  accuracy.  Having  a  baseline  test  will  eventually  allow  the  BRL- 
CAD  ray-tracer  to  yield  comparative  results  that  could  be  assessed  relative  to  itself  or  to  other 
ray-tracers  (such  as  FASTGEN,  a  ray-tracer  popular  in  many  US  Air  Force  applications). 


2.  Experiment/Calculations 


This  project  was  approached  by  modifying  2  example  programs  in  the  BRL-CAD  source  code. 
The  first  program,  wdb  example.c,  creates  a  database  of  BoTs  by  reading  information  about  the 
triangles  from  a  file,  wdb  example.c  was  modified  to  create  the  database  of  BoTs  for  this 
project.  The  second  program,  rtexample.c,  implements  a  sample  ray -tracer  by  loading  a  triangle 
from  a  database  and  firing  a  single  ray  toward  a  point  relative  to  the  triangle  and  in  a  direction 
normal  to  the  triangle.  The  rtexample.c  was  modified  to  fire  rays  at  a  multitude  of  points  relative 
to  each  triangle  in  the  database  and  at  various  tolerance  levels. 

The  wdb_example.c  was  used  as  a  starting  point  and  modified  to  create  a  new  program, 

(wdb  triangle.c)  that  would  generate  all  possible  BoTs.  This  process  required  several  steps. 
Initially,  wdb  triangle.c  was  modified  to  create  a  database  of  a  single  triangle  using  user 
specified  coordinates  of  the  3  vertices  input  from  the  command  line.  This  program  was  then 
expanded  to  create  all  possible  triangles  for  a  bounding  box  and  a  step  size  in  the  x,  y,  and  z 
directions  specified  by  the  user.  The  BRL-CAD  library  function  for  creating  BoTs,  mk  bot,  was 
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continually  called  during  each  iteration  of  the  9  nested  for-loops  in  the  main  function  to  create  a 
triangular  object  in  the  database  for  each  possible  combination  of  vertices.  This  resulted  in 
multiple  objects  in  the  database  that  had  the  appearance  shown  in  Fig.  1  when  opened  using  the 
graphics  interface  application,  MGED. 


0  0  0  \  MCED  7.25.0  Graphics  Window  (id_0)  - /Users/Charith_Ranawake/brlcad/build/tri.g  -  Upper  Ri... 

File  Edit  Create  View  ViewRing  Settings  Modes  Mjsc  Tools  Help 


Fig.  1  All  possible  triangles  within  bounding  box  (1,  1,  1)  with  step  size  of  0.5 

A  ray-tracing  program,  rt  bot.c,  was  then  implemented  to  fire  rays  at  several  locations,  both  on 
and  off  the  triangle,  for  every  triangular  object  in  the  database  that  was  created  by 
wdb  triangle.c.  These  locations  are  shown  in  Fig.  2. 
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Fig.  2  Rays  fired  at  triangle 

Figure  3  shows  the  basic  outline  of  the  ray-tracing  program,  rtbot.c.  The  main  program  has  the 
same  basic  structure  as  the  corresponding  wdbtriangle.c  program  that  created  the  BoTs.  The 
user  inputs  the  bounding  box  and  the  step  size  of  the  program  as  runtime  arguments.  The 
program  uses  9  nested  for-loops  for  incrementing  vertex  coordinates.  For  each  combination  of 
vertices  that  forms  a  triangle,  wdb  triangle.c  loads  the  corresponding  object  from  the  BoT 
database  and  calls  a  function  for  performing  the  ray -tracing  on  the  triangle.  As  mentioned  before, 
ray-tracing  involves  calculating  the  coordinates  of  the  target  point,  the  unit  normal  to  the 
triangle,  and  the  coordinates  of  the  firing  point.  The  ray-tracing  function  also  maintains  counters 
to  hold  such  information  as  the  number  of  rays  fired  and  the  number  of  successes  (hits  or 
misses).  For  interior  target  points,  vertices,  and  edges,  a  success  occurs  when  the  ray-tracer 
records  a  hit.  For  an  exterior  point,  a  success  occurs  when  the  ray-tracer  records  a  miss. 
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Initialize  variables  for  ray-tracer 


No 


> 


Set  vertices  of  a  triangle 
Load  object  from  database 


> 


Calculate  target  point 
Calculate  direction  and  origin  of  ray 
Fire  ray  and  increment  count  if  success 


Fig.  3  Ray-tracing  program  flowchart 


Several  steps  are  involved  when  determining  the  x,  y,  and  z  coordinates  for  a  target  point.  Figure 
4  shows  an  example  in  which  the  target  point  is  interior  to  the  triangle  and  on  the  median  joining 
vertex  B  to  mid-point  of  line  AC.  The  first  step  is  to  calculate  the  midpoint,  D  of  line  BC,  as  the 
target-point  coordinates: 


(xi+x3)  _(yx+y3)  „  _(zx+z3) 


->y* 


,Z4 


(1) 


Then  the  distance  BD,  which  we  will  call  “L”  the  distance  to  the  target  point,  is  calculated  from 
using  the  formula  below: 


L  =  J(x2-xJ+(y2-yJ+(z2-zJ 


(2) 
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Fig.  4  Target-point  diagram 

The  target  point  is  calculated  using  Fig.  5  by  applying  the  following  proportion  formula: 


m 

JL 


O 


nr 

-o- 


n 

JL 


1 

O 


B  (x2,  y2,  z2) 


D  (x4)  y4,  z4) 


Target  (xt,  yt,  zt) 


Fig.  5  Coordinates  of  a  point  dividing  a  line  segment 

Where  the  target-point  coordinates  (xt,  yt,  zt)  are  calculated  by  Eq.  3  using  the  proportional 
distances  m  and  n  relative  to  L,  where  m-L-n,  and  n  =  tolerance 

mx,  +  nx7 

xt=— - 1 

m  +  n 
my4  +  ny2 

yt= — - —  .  C 

m  +  n 
mz,  +  nzn 

zt=~ - ^ 

m  +  n 

These  calculations  are  specifically  for  points  that  are  interior  to  the  triangle.  The  exterior  points 
can  be  calculated  in  a  similar  manner  using  a  slightly  different  formula. 
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In  addition  to  the  target  point,  the  firing  direction  and  the  firing  point  were  calculated.  Figure  6 
shows  the  case  where  a  ray  is  fired  from  the  firing  point  (xf,yf,zf)  in  a  direction  normal  to  the 

triangle. 


(xf  yf ,  zf) 


C(x3,  y3,  z3) 


Fig.  6  Triangle  with  normal  vector 


Using  the  coordinates  of  the  target  point,  the  firing  point  and  the  direction  of  the  ray  were 
calculated  as  follows: 

•  Step  1 :  Calculate  normal  to  AABC  using  the  cross  product  of  2  vectors  (Ex:  AB  x  AC). 

•  Step  2:  Calculate  the  unit  vector  n  =<  nx,ny,nz  >of  the  cross  product. 

•  Step  3:  Determine  the  coordinates  of  the  firing  point  (xf,y  zf)  as  the  following: 


xf  =  x,  —  hnx 
yf=y,-hny 
Zf=zt-  hnz 


(4) 


3.  Results  and  Discussion 


To  test  the  performance  of  the  ray-tracer,  a  BoT  was  created  for  the  bounding  box  (2,  2,  2)  with 
the  step  size  of  1.5.  This  created  a  database  of  336  triangles.  Then  the  ray-tracer  program  was 
executed  on  the  set  of  triangular  objects;  the  number  of  successes  and  the  number  of  attempts 
were  recorded  for  each  type  of  target  point.  The  data  are  tabulated  in  Table  1 .  Figure  7  is  a  bar 
chart  showing  the  variation  of  percentage  of  successes  with  the  type  of  target  point.  The  data 
show  that  the  ray-tracer  recorded  100%  success  for  interior  points  as  expected.  Also,  as  the 
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tolerance  increased  from  0.00001  to  0.1,  the  percentage  of  successes  increased  from  91%  to 
95%.  The  ray-tracer  recorded  the  smallest  percentage  of  successes  when  the  target  point  was  on  a 
vertex  or  edge.  From  these  results,  it  can  be  concluded  that  BRL-CAD’s  ray-tracer  is  fairly 
accurate  for  target  points  that  are  interior  to  the  triangle  as  well  as  target  points  that  are  small 
tolerance  levels  away  from  the  edges  and  vertices.  However,  the  accuracy  decreased  as  the  rays 
were  fired  at  a  target  point  that  was  on  a  vertex  or  edge. 

Table  Data  for  the  ray-trace  of  BoTs  for  the  bounding  box  (2,  2,  2)  with  step  size  1.5 


Database:  Trian 

gles  (336) 

Type 

Successes 

Attempts 

Percentage  (%) 

Interior 

1344 

1344 

100 

Vertex/Edge 

1710 

2016 

84.82142857 

Tolerance  -  0.1 

5802 

6048 

95.93253968 

Tolerance  -  0.001 

5772 

6048 

95.43650794 

Tolerance  -  0.00001 

5528 

6048 

91.4021164 

Database:  TRIANGLE  (336  Triangles) 


100.0000 

90.0000 

g  80.0000 

S  70.0000 
1/1 

d  60.0000 
£  50.0000 

O 

m  40.0000 

ro 

U  30.0000 

u 

£  20.0000 
10.0000 
0.0000 


Interior  Vertex/Edge  Tolerance  -  0.1  Tolerance  -  0.001  Tolerance - 

0.00001 


Location 


Fig.  7  Data  for  the  percentage  of  successes  for  the  bounding  box  (2,  2,  2)  with  step  size  1.5 
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4.  Summary 


In  this  project,  BRL-CAD’s  ray-tracer  was  tested  on  the  primitive  Bag  of  Triangles  (BoTs).  Two 
programs  were  created  to  accomplish  these  tests.  The  first  program  created  numerous  triangles  or 
objects  in  the  database,  while  the  second  program  was  used  to  fire  rays  at  each  object.  Based  on 
the  results,  the  ray-tracer  was  very  accurate  on  points  on  the  interior  of  the  triangles  and  small 
tolerances  away  from  the  vertices  and  edges.  But  the  accuracy  of  the  ray-tracer  decreased 
significantly  when  tested  on  edge  or  vertex  points. 
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